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Abstract 

We study the physical propagating modes in a massive gravity model in curved cos- 
mological backgrounds, which we have found as classical solutions in our previous paper. 
We show that, generically, there exist such the cosmological background solutions con- 
sistent with the equations of motion where we assume the ghost condensation ansatzes. 
Using the (l+3)-parametrization of the metric fluctuations with 'unitary' gauge, we 
find that there is neither a scalar ghost nor a tachyon in the spectrum of the propagat- 
ing modes, the tensor modes become massive owing to gravitational Higgs mechanism, 
and the model is free of the Boulware-Deser instability. The price we have to pay is 
that the scalar sector breaks the Lorentz-invariance, but there are no pathologies in the 
spectrum and lead to interesting phenomenology. Moreover, we present a proof of the 
absence of non-unitary modes for a specific ghost condensation model in a cosmological 
background. 
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1 Introduction 



In recent years, we have watched revival of interests to construct massive gravity theories from 
different physical motivations [1]-[16]. It concerns a very simple question: Can the graviton 
have a (small) mass consistent with experiments? If so, how? 

One motivation behind this question comes from the astonishing observational fact that 
our universe is not just expanding but is at present in an epoch of undergoing an accelerat- 
ing expansion [17, 18]. Although the standard model of cosmology is remarkably successful 
in accounting for many of observational facts of the universe, it is fair to say that we are 
still lacking a fundamental understanding of the late-time cosmic acceleration in addition to 
problems associated with dark matter and dark energy. 

Several ideas have been thus far put forward for explaining this perplexing fact [19]. It is 
here that massive gravity theories might play an important role since massive gravity theories 
could modify Einstein's theory of general relativity at large cosmological scales and might lead 
to the present accelerated expansion of the universe without introducing still mysterious dark 
matter and dark energy at all. Since general relativity is almost the unique theory of massless 
spin 2 gravitational field whose universality class is determined by local symmetries under 
general coordinate transformations, any infrared modification of general relativity cannot 
aviod introduction of some kind of mass for the graviton. 

The other motivation for attempting to construct massive gravity theories lies in string 
theory approach to quantum chromodynamics (QCD) [8]. For instance, as inspired in the 
large-N expansion of the gauge theory, which defines the planar (or genus zero) diagram 
and is analogous to the tree diagram of string theory, if we wish to apply a bosonic string 
theory to the gluonic sector in QCD, massless fields such as spin 2 graviton in string theory, 
must become massive or be removed somehow by an ingenious dynamical mechanism since 
such the massless fields do not appear in QCD. Note that this motivation is relevant to the 
modification of general relativity in the short distance region while the previous cosmological 
one is to the infrared modification. It is worthwhile to point out that the modification of both 
short and long distance regions reminds us of T-duality R <-> l 2 s / R in string theory where l s 
is the fundamental string length scale. If there is the interesting symmetry behind the both 
distance regions in general relativity, such a modification could perhaps open a fruitful dialog 
between cosmology and QCD. 

An approach for the construction of massive gravity theories is to take account of the 
spontaneous symmetry breakdown (SSB) of general coordinate reparametrization invariance 
[l]-[8]. Note that the SSB of general coordinate reparametrization invariance might provide a 
resolution for cosmological constant problem [6, 8]. Though the cosmological constant problem 
needs a theory of quantum gravity, solving the problem requires a low energy mechanism, for 
instance, a partial cancellation of vacuum energy density stemming from quantum loops. 
We expect that from the analogy with the Higgs mechanism in conventional gauge theories, 
the Nambu-Goldstone mode mixes with the massless graviton, thereby changing the vacuum 
structure of gravitational sector in a non-trivial way. 

Recently, as a device for performing a consistent infrared modification of general relativity 
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an idea of ghost condensation has been proposed [9]. In this scenario, the 'unitary 1 propagat- 
ing scalar field appears as the Nambu-Goldstone boson for a spontaneously broken time-like 
diffeomorphism, and yields a possibility for resolving cosmological problems such as inflation, 
dark matter and dark energy. The key observation here is that a scalar ghost is converted to 
a normal scalar with positive-energy excitations. 

More recently, 't Hooft has proposed a new Higgs mechanism for gravity where the mass- 
less graviton 'eats' four real scalar fields and consequently becomes massive [8]. In his model, 
vacuum expectation values (VEV's) of the scalar fields are taken to be the four space-time 
coordinates by gauge-fixing diffeomorphisms, so the whole diffeomorphisms are broken spon- 
tanously. Of course, the number of dynamical degrees of freedom is left unchanged before 
and after the SSB. Actually, before the SSB of diffeomorphisms there are massless gravitons 
of two dynamical degrees of freedom and four real scalar fields whereas after the SSB we have 
massive gravitons of five dynamical degrees of freedom and one real scalar field. Afterward, 
a topological term was included to the 't Hooft model where an 'alternative' metric tensor is 
naturally derived and the topological meaning of the gauge conditions was clarified [20]. 3 

The problem in the 't Hooft model is that a scalar field appearing after the SSB is a non- 
unitary propagating field so that in order to avoid violation of unitarity it must be removed 
from the physical Hibert space in terms of some procedure. 4 A resolution for this problem 
was offered where one requires the energy- momentum tensor of the matter field to couple to 
not the usual metric tensor but the modified metric one in such a way that the non-unitary 
scalar field does not couple to the energy-momentum tensor directly [8]. 

In this context, one could imagine that the ghost condensation scenario enables us to avoid 
emergence of the non-unitary scalar field in gravitational Higgs mechanism since the annoying 
non- unitary scalar field has its roots in one time-like component (i.e., ghost) in four scalar 
fields, to which the ghost condensation idea could be applied. Thus, it should be regarded 
that this observation is an alternative method for removing the non-unitary propagating scalar 
field in gravitational Higgs mechanism by 't Hooft. 

In our previous article [11], as a first step for proving this conjecture, we have studied 
classical solutions in the unitary gauge in general ghost condensation models. This analysis is 
needed for understanding in what background gravitational Higgs mechanism arises. It turns 
out that depending on the form of scalar fields in an action, there are three kinds of classical, 
exact solutions, which are (anti-) de Sitter space-time, polynomially expanding universes and 
flat Minkowski space-time. 

In this article, we wish to prove that there is indeed no non-unitary mode in the spectrum 
of propagating modes around cosmological backgrounds and the graviton becomes massive 
because of gravitational Higgs mechanism in these models. Thus, the models at hand are free 
of the problem associated with the non-unitary propagating mode and the Boulware-Deser 
instability [22] in curved backgrounds. 

3 Similar but different approaches have been already taken into consideration in Ref. [21]. 

4 More recently, models of gravitational Higgs mechanism without the non-unitary propagating scalar field 

were proposed in Ref. [10]. 
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This paper is organized as follows: In section 2, we review that there is an expanding 
cosmological solution with zero acceleration to the equations of motion in a massive gravity 
model, what we call, 't Hooft model with ghost condensation. In section 3, we consider a 
more general model and look for classical solutions where we will find expanding cosmological 
solutions with non-zero acceleration. In section 4, based on the massive gravity model in 
section 3, the propagating modes around the cosmological backgrounds are examined in detail 
by using the (l+3)-parametrization of the metric fluctuations. We find that the tensor modes 
become massive, the vector modes are non-dynamical, three of the scalar modes are also 
non-dynamical and one scalar becomes massive. This scalar is originally non-unitary mode, 
that is, a ghost, but becomes a normal particle because of the ghost condensation mechanism. 
However, the dispersion relation is not usual one, so this mode breaks the Lorentz invariance 
as expected from the ghost condensation scenario. In section 5, we present a proof of the 
absence of non-unitary modes for the specific ghost condensation model treated in section 2. 
The final section is devoted to conclusions and discussion. 



2 Cosmological solution in 't Hooft model with ghost 
condensation 

In this section, we wish to review our previous work [11] and explain how to derive classical 
solutions in massive gravity models. 

Let us start with review of our 't Hooft model with ghost condensation. The general proof 
that there is no non-unitary scalar ghost in this model will be given in section 5. The action 
of ghost condensation [9], which is inspired by gravitational Higgs mechansim, takes the form 

5. 

S = -^-Jd A x^[R-2K + f{X)-g^d,rdA b ria^ 

= j d A x^[R - 2A + F(X) - g^drfdrfSij], (1) 

where r] a b is the internal flat metric with diagonal elements (— 1, +1, +1, +1). The indices 
/i, v[= 0,1,2,3) and 1,2,3) are space-time and spatial indices, respectively. And we 

have defined F(X) = X + f{X) where f(X) is a function of X which is defined as 

X = g^d^d^ . (2) 

Let us note that the last term in the first equality in (1) is added to the original ghost 
condensation action, or alternatively speaking from the side of gravitational Higgs mechanism, 
the third term f(X) is introduced to the original action of gravitational Higgs mechanism. 

5 For simplicity, we have put the Newton constant Gat at the front, so the dimension of ghost-like scalar 
field <jP differs from that of the original ghost condensation model, but it is easy to modify the dimension by 
field redefinition. 
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Recall that the ghost condensation scenario consists of three ansatzes: The first ansatz 
amounts to requiring that the function F{X) has a minimum at some point X min such that 

Fx(X min ) = 0, Fx X {X m in) > 0, (3) 

where Fx{X min ), for instance, means the differentiation of F{X) with respect to X and then 
putting X = X min . Note that the latter condition in (3) ensures ghost condensation. 

As the second ansatz, the ghost-like scalar field 0° is expanded around the background 
mt as 

0° = mt + 7T, (4) 

where 7r is the small fluctuation. This equation can be interpreted as follows: Using a time- 
like diffeomorphism <50° = e°, one can take the gauge tt = 0. In other words, n is a Nambu- 
Goldstone boson associated with spontaneous symmetry breakdown of time traslation. 

The last ansatz is natural from the second ansatz. That is, since time coodinate t plays a 
distinct role from spatial coordinates x l {i = 1, 2, 3), it is plausible to make an assumption on 
background space-time metric, which is of the Friedmann-Robertson- Walker form 

ds 2 = -dt 2 + a{t) 2 dQ 2 , (5) 

where a{t) is the scale factor and dfl 2 is the spatial metric for a maximally symmetric three- 
dimensional space. 

The equations of motion are easily derived as follows: 

d^y/^gg^drf) = 0, 
d,{^g^Fx{X)d v 4> ) = 0, 

G^ + Ag^ = T^, (6) 

where the stress-energy tensor is defined as 

= {d^dA* ~ \g^d^ a d^ h )r]ai 

- {d,<tPd v <tPf x {X)-\g, v f{X)). (7) 

Now let us solve these equations of motion under the ansatzes of ghost condensation and 
the unitary gauge for diffeomorphisms 

a = mx^d^. (8) 

It is easy to check that the 0*-equations are trivially satisfied. Moreover, the 0°-equation 
reduces to the expression 

d Q {a(tfF x {X)) = 0, (9) 
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whose validity requires us that one of ghost condensation ansatzes should be automatically 
satisfied, i.e., F x (X min ) = since X min = X(0° = mt) = —m 2 and d X min = 0. 
Finally, Einstein's equations are cast to 

,d, 9 ~ 3m 2 
a ~ = 2tf; 

_ 2 £ _(£)* + A = (10) 



a 



where we have defined A = A — \F{ — m 2 ) and the dot over a(t) denotes the derivative with 
respect to t. Here we have made use of the expression of the stress-energy tensor which is 
obtained by inserting Eq. (8) to Eq. (7) 

Tfiv = m\^-\g^9 ab Va b )+m 2 5l^ + \g»vf{-m 2 )- (11) 
From the equations (10), we can eliminate the terms involving a whose result is written as 

3- = A. (12) 
a 

It then turns out that the general solution for this equation exists only when A = 0. Thus, 
Eq. (12) reads a = 0, so we have the general solution 

a(t) =c(t-t ), (13) 

where c and t are integration constants. Substituting this solution into the first equation in 
Eq. (10), c is fixed to c = ±^|. As a result, the line element takes the form 

ds 2 = -dt 2 + ^-(t-t ) 2 dn 2 , (14) 

which describes the linearly expanding universe with zero acceleration. 6 Notice that this 
solution is a unique classical solution, so a flat Minkowski space-time, for instance, is not a 
solution of this model. Here it is worthwhile to mention that although the present universe 
seems to be accelerating so that the linearly expanding universe is excluded from WMAP 
experiment [23], the general solution (14) might be useful in describing the future or past 
status of universe. 



3 Cosmological solutions in more general massive grav- 
ity models 

In the previous section, we have found an interesting cosmological solution, which is very 
similar to the Milne universe, in the 't Hooft model with ghost condensation. However, the 

6 This solution was also obtained in a different massive gravity model [10]. 
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solution has zero acceleration so that it does not describe the present epoch of the acceler- 
ating universe. It would be more desirable if we could get a classical solution with non-zero 
acceleration. It has been shown in our previous paper [11] that this desire is actually real- 
ized if one generalizes the model in section 2 to, what is called, a general ghost condensation 
model. In this article, instead of explaining the derivation of the solutions in the general ghost 
condensation model, we show that the solutions also exist in more general massive gravity 
models with the ghost condensation potential. 

In this section, we will obey the following line of arguments: In order to make the analysis 
of propagating modes easier, which will be done in the next section, we will first choose a 
simple massive gravity model in a general class of those models. Next, using this specific 
model, we will show that there are polynomially expanding universes' solutions with non-zero 
acceleration. Finally, we will examine under what conditions the solutions at hand remain 
the solutions to the equations of motion in the most general models. 

Now we shall take a different line element from (5) since the form is more convenient for 
later analysis in section 4. Of course, at a final stage it is easy to transform to the expression 
like (5) through a global coordinate transformation. The form we take as the line element is 
of conformal flat type: 

ds 2 = g liV dx il dx u 

= e 2A{t) r]^dx^dx u , (15) 

where rj^ is the flat metric. 

We consider a simple massive gravity model 7 : 

S = - [ ^-Jd 4 x v ^[R-2A + F(X)-G(W^)} 

= Ta~7^~ I d 4 x[L E H + L A + L F + L G ], (16) 

where X is defined in (2) and W tj is defined via new variables Y' lj and V 1 as 

Y'YJ 

w»j = Y l ° - — — 
X ' 

Y ij = gTd^W, 

V 1 = gTd^drf. (17) 
This model of massive gravity was obtained by considering the residual diffeomorphisms 

x* -» x H = x i + C(t), (18) 



7 For the analogy with the previous model in section 2, we have introduced two different functions F(X) 
and G(W y ') separately for variables X and W y ', but it is more economical to do only one function F(X, W l i) 
dependent on two variables. This simple model has been investigated in Refs. [14]-[16]. 
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which are translated to the symmetries of the scalar fields as 

0*-0 ,i = 0*+aA (19) 

where 0(t) and £ l (0°) are the infinitesimal transformation parameters which are dependent 
on t and 0°, respectively. Note that both X and are invariant under (19). 

Let us show that there are polynomially expanding cosmological solutions with non-zero 
acceleration in this model. To this end, we first derive the equations of motion. The 0*-, 0°-, 
and Einstein's equations are derived in a straightforward manner 

yi yiyi fid 

d»[V^9^m F x + ( Y d„<f> i -—d„ ( f ) )—}} = 0, 

G^ + Ag^ = T M „, (20) 

with the stress-energy tensor being given by 

2 yiyj Be 

+ (21) 
With the unitary gauge (8), each variable takes the form 

X = m V° = -mV 2i = l, 

V* = m 2 g 0i = o = v i , 

= m 2 g 13 = m 2 e- 2A Sij = Y lj , 

W ij = m 2 e~ 2A ^ = WSij = W ij . (22) 

Then, the ^-equations are reduced to 

^«H0, (23) 

which is trivially satisfied since the quantity in the square bracket depends on only t. 
Next, the 0°-equation reads 

d t [e 2A F x (X)} = 0, (24) 

which is satisfied by one of ghost condensation ansatzes, F X (X) = 0. Furthermore, Einstein's 
equations are recast to 

A(-25y u - 2^) + (A) 2 (25y u - Vfiu ) + Ae 2 V 
= m X 6 l^^ + r 2A ^AF(X) - G(W^)}. (25) 
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In order to proceed further, we need to fix the form of the potential term G(W 1 ^). It then 
turns out that there are non-trivial cosmological solutions provided that we select 

G(W ij ) = l -KTr(W li ) n - 2A, (26) 
o 

where K is a positive constant and A = A — ^F(X). For later convenience, using this 
expression (26) let us define the following two quantities 

aw* v 

d 2 G(W l i) 
dW^dW kl 

where G\ and G2 are defined as 



G 2 (SikSji + SuSjk), (27) 



G l = ^nW n - l = ^n{m 2 e~ 2A ) n -\ 

G 2 = ^^^W^ = ^^^(rn 2 e- 2A r 2 . (28) 

O Z O Z 



As a result, Einstein's equations read 



(A) 2 = -l[F(X) -G(W ij ) -2A]e 2A , 
6 

A+\(A) 2 = -UF(X)-G(W ij )-2A]e 2A -lm 2 G 1 . (29) 
These equations are easily integrated to be 

e A = [±M(n-l)m n (t-t )]^, (30) 
V 6 

where t is an integration constant. Then, the line element becomes (after a suitable redefi- 
nition of by an overall constant factor) 

ds 2 = (t - t Q )^r] ilv dx tJ, dx v 

= —dr 2 + (r — To)"dx t dx i , (31) 

where the latter expression informs us that the second derivative of the scale factor, that is, 
the acceleration, is positive for n < 1. Of course, this condition might not be meaningful 
phenomenologically since the present universe seems to be entering a new era of exponentially 
accelerating inflation again, which is controled by the equation of state w — ^ — — 1, whereas 
our solutions describe the polynomially accelerating universes. Nevertheless, we shall take 
account of this condition since at present there do not exist sufficient evidences to exclude 
the polynomially accelerating universe such as quintessence (—1 < w < — |) and phantom 
(w < —1) etc. 
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In this way, we have obtained an interesting class of cosmological solutions, which have 
a behavior of polynomially expanding universes with non-zero acceleration, by generalizing a 
massive gravity model in such a way that the potential terms in the action include functions 
of not only 0° but also 0\ 

To close this section, it is valuable to ask ourselves whether such the solutions exist even 
in the most general models or not. Recall that the symmetries (19) constrain the model to 
some degree in the sense that the potential terms are not a general function of X, V 1 and Y 1 ^ 
but the more restricted function of X and which are invariant under (19). 

We therefore take the potential term to be the most general function, which is an arbitrary 
function of X, V i and Y ij : 

S = J d'x^-g[R - 2A + U(X, V\ Y ij )]. (32) 

Now let us derive the equations of motion under the unitary gauge (8) whose concrete ex- 
pressions are given by 

d [^(\g^^L + griJ^L)] = 0, 

0^9(9^ + \r§;)] = 0, 

G^ + Ag^ = T^. (33) 
Now, the second 0°-equation is satisfied if one imposes the ghost condensation ansatz 

ajHx^l = o. (34) 

Next, we find that the first 0*-equations are satisfied if and only if there exist new conditions 
on U 

dU{X,V\Y ij ) _ 

In this regard, note that in order to keep the 5*0(3) rotational symmetry, U must be a 
function of V l V % , V l YijV 3 etc., so 9U gyi ^ would be proportional to V 1 . In the unitary gauge, 
V 1 is vanishing so that the conditions (35) are valid. 

The remaining equations of motion are Einstein's equations, which read in the unitary 
gauge 

A(-2$fi - 2^) + (A) 2 (26°X - n,,,) + Ae 2 V 

= - m>W^ + Sl6l§- + W^L). (36) 

The requirement that these equations should have the same solutions as (30) amounts to the 
similar conditions to (26) and (28), which are 

U(X, V\ Y ij ) = -K(m 2 e- 2A ) n + 2A, 
dU(X,V\Yv) _ * , ,-2Ax„-i 



dYv 3 



_ n(m V^)»-M y . (37) 
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These conditions appear in a natural way when we assume the potential term U to be the 
polynomial type like (26), so we could conclude that the classical solutions (30) exist as well 
in the most general massive gravity models by picking up an appropriate form of the ghost 
potential. 



4 Analysis of propagating modes 

In this section, we analyse the physical modes propagating in the cosmological backgrounds 
obtained in section 3. Thus, let us consider small fluctuations around the metric 

g l u, = e 2A( - t \r tl w + h l »,). (38) 

With Eq. (38), some quantities relevant to the metric tensor read in the second-order ap- 
proximation level of the fluctuations 



1 1 

2 4 ^ ' 8 



yfZjjgi" = e 2A{t) [rf - hT + \rf v h + {-]h aP h a(3 + -h 2 ) V ^ 

2 4 8 

- -hh^ + h^ a h u a ], (39) 

etc. Here, summation over space-time indices \i = 0,1,2,3 is carried out with the flat 
Minkowski metric rj^ while that over spatial indices i — 1,2,3 is done with the Kronecker 
5ij, so we have, for instance, 

h = = -/loo + K. (40) 

In the analysis of the physical modes, it is convenient to make use of the (1 + 3)- 
parametrization of the metric fluctuations [24] 

hoa = 20, 

hoi = Si + diB, 

hij = hi? -(d i F j + d j F i )-2(i>5 ij -d i d j E), (41) 



where hj? is a traceless, transverse spatial tensor 



dih^ = = h™, (42) 
and Si and Fj are transverse spatial vectors 

diSi = = d t F t . (43) 
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In this parametrization, the tensor hf? has 2 degrees of freedom (d.o.f.), the vectors Si and 
Fi do 2 x 2 = 4 d.o.f., and the scalars 0, 5, ^ and E do 1 x 4 = 4 d.o.f., so in total 10 d.o.f. 
which is exactly equal to the number of h^. 

Incidentally, recall that in Einstein's general relativity only the massless spin 2 graviton 
hj? of 2 d.o.f. is physical owing to general coordinate invariance. On the other hand, in a 
general massive gravity model 8 , the spin 2 tensor hj? of 2 d.o.f., the spin 1 vector Fi of 2 d.o.f. 
and the spin 1 scalar ip (or E) of 1 d.o.f. are physical so totally we have the massive graviton 
of 5 d.o.f. (the other modes are non-dynamical) which involves spins ±2, ±1 and 0. The 
important point to notice is that there in general remains one scalar mode E (or ip), which 
has negative norm, so it is in essence a ghostl This peculiar feature, that is, the existence of a 
scalar ghost in the physical Hilbert space, is usually avoided by removing the ghost with the 
help of an enhancement of gauge symmetry as in the Fierz-Pauli Lorentz-covariant massive 
gravity model [25]. However, this nice property is lost in curved backgrounds owing to the 
disappearance of the enhanced gauge symmetry, and in consequence we have the Boulware- 
Deser instability scalar mode [22]. 

Accordingly, the real problem is to find a method of removal of the ghost mode from the 
physical spectrum without recourse to the enhanced gauge symmetry. In this article, we try 
to remove the ghost by appealing to the ghost condensation mechanism. 

After a little lenghty calculation, the quadratic part of the Einstein-Hilbert action and 
the cosmological term in (16) is decomposed into tensor, vector and scalar sectors up to total 
derivative terms 

T (2) _ T (T) , T (V) T (S) 

^EH — ^EH "r ^EH "r ^EHi 

Lf = 4 T, + 4"> + £!?, (44) 



where each term is given by 

4 T i = e^[-{i+i(i) 2 }(0 2 +Umn 2 -mrf}i 



L%1 = e 2A [3(A) 2 S 2 + 2{A + ^(A) 2 }F t AE l --(S t + d F i )A(S t + d F t )}, 



1 •. 1 



Leh = e 2A [(A) i {-90 2 -9^ + 2(90^-30AE + ^(^ J B) 2 )-2(M^+ )} 
+ A(A<pAB + 120<9 O ^ - 4(j)d AE) - 4i(^A£ + -(AE) 2 ) 

- 4doi/>AB - 4<i>Ail> + 6il>dfy -24>Ai> + 4d il>d AE], (45) 



and 



4 T) = \^\h^)\ 

L { P = -Ae 4A (^ + F J AF J ), 

Lf = -2Ae^[-i0 2 + 30^-0AE + ^ 2 -M^-^(A^) 2 + i(^) 2 ], (46) 



3 Here the term general means that there is no special residual gauge symmetry. 
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where A denotes the Laplacian operator defined as A = df. 

Furthermore, provided that by L F+G we denote the sum of the quadratic part of the 
potential terms Lp and Lq, it is also decomposed into each spin sector as 

L?l G = [-\e 4A {F-G)-e 2A m 2 G 1 -m 4 G 2 ]{h^) 2 , 

4+g = [\e 4A (F-G) + 2e 2A m 2 G 1 + 2m 4 G 2 ]F l AF t + ^e 4A (F-G)S 2 , 

Lf +G = e 4A [-±(F-G) + 2F xx X 2 }<f> 2 + ±e 4A (F-G)(d i B) 2 

+ [e 4A (F -G) + 2e 2 V 2 G 1 ]0(3^ - AE) + [^e 4A (F - G) + §e 2A m 2 G x - \2m 4 G 2 \i\) 2 
+ \—e 4A [F - G) - Ae 2A m 2 G 1 + 8m 4 G 2 }^AE 

+ [-- e 4A (F-G)-2e 2A m 2 G l -4m 4 G 2 ](AE) 2 , (47) 
2 

where F, F xx and G denote F(X), F XX (X) and G{W ij ), respectively. 

We are now in a position to discuss each spin sector in order. In the tensor sector, the 
total Lagrangian takes the form 

r(T) _ AT) f (T) AT) 
^ = ^EH + -^A + ^F+G 

= \e 2A [(d hJ?) 2 - {d^f] + [-\e 2A m 2 Gl - m 4 G 2 }(^) 2 , (48) 

where the background equations, those are, Einstein's equations (29) are used to simplify the 
expression. Then, the equations of motion for hj^ give 

nhl? - 2Ad hJ? - 2(m 2 Gi + 2m 4 G 2 e- 2A )k[T = o, (49) 

where the d'Alembertian operator is defined as □ = <9 M <9 M = — d 2 + A. 
In order to see what mass of the graviton is, let us introduce 

h = e A ^J. (50) 

In terms of hij, the equations of motion (49) read 

Dhij + [A+ (A) 2 - 2{m 2 G 1 + 2m 4 G 2 e~ 2A )\h lj = 0. (51) 

Thus, the effective mass square of the tensor modes, M 2 is given by 

M 2 = -[A + (A) 2 -2(m 2 G 1 + 2m 4 G 2 e~ 2A )} 
An 2 + n - 2 1 
(n-iy (t-t ) 2 ' 

where we have used (30) and (28) in the second equality. 



(52) 
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Thus, the tensor modes in the cosmological backgrounds are massive for n > or 
n < — (n 7^ 1) while for ^g 7 ^ < n < they are effectively tachyonic and the 

backgrounds become unstable at least perturbatively. For n = , the graviton becomes 

massless. Moreover, the effective mass of the graviton approaches zero in the limit t — > oo. 

We next move to the vector sector. With the help of Einstein's equations (29), the total 
Lagrangian for vector modes can be written as 

= + + L F + G 

-- "^(3 + d F t )A(S t + doF,) + {e 2A m 2 G l + 2m 4 G 2 )F l AF l . (53) 
Taking the variation with respect to Si, we have the equation 

Si = -d F % . (54) 
Substituting it into L^ v \ we have only the second term in (53) 

LW = {e 2A m 2 G 1 + 2m 4 G 2 )F i AF l . (55) 
The equations of motion for F$ therefore yield 

Fi = 0. (56) 

Plugging it back into in (55) again, the Lagrangian becomes identically vanishing, so 
that Fi do not obey any equations of motion and take any values. 

This arbitrariness, of course, is a consequence of the residual diffeomorphism invariance 
(18). In this case, the reparametrization symmetries read 

5 V = 2i^C° + d„Cu + 0„Ci. (57) 

With d = d(t) in (18) and the parametrization (41), we have the residual gauge symmetries 

5S t = d d, SF t = -d, (58) 

from which the modes Fi and Si become non-dynamical. More precisely speaking, it is Fi 
that the residual gauge symmetries (58) make non-dynamical since Si remain non-dynamical 
irrespective of the existence of mass term as easily seen in (53). In this way, we have proved 
that all the vector modes are not physical but non-dynamical in this simple massive gravity 
model. 

Finally, we are ready to examine the scalar sector, which is known to be the most prob- 
lematic and complicated. In fact, the Boulware-Deser mode [22] appears in this sector in the 
Fierz-Pauli massive gravity model [25]. 
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After utilizing Einstein's equations (29), the total Lagrangian is of form 



r (S) _ AS) j (5) AS) 

^ = ^EH + -^A + ^F+G 

= 2e 2A [-20A^ - 2d ipAB + 2d ipd AE + S^dfy - ipAip 

+ A(2(f)AB - 2(f)d AE + 60<9 O VO] + e 4A (F -G-2A + 2F xx X 2 )<f) 2 

+ 3[e 4A (F - G - 2A) + 2e 2A m 2 G 1 - 4m 4 G 2 ]^ 2 + 2m 2 G 1 e 2A (30^ - 0A£) 

- 2m 2 G'ie 2A [^A^ + ^(A^) 2 ] + 8m 4 G' 2 [M^-^(A J B) 2 ]. (59) 
The equation of motion for B gives 

d ip = A<f). (60) 

Thus, integrating over B and <fi and using Einstein's equations again, up to total surface terms 
j s reduced to 

= -^—m 4 F xx d ^d ^ + 2ne 2A ^AiP-3nm 2 e 2A G 1 ^ 2 
(A) 2 

_ le^d^AE - 2e 2A m 2 G 1 [^AE + -(AE) 2 ] 

+ 8m 4 G 2 [V'Ai? - i(AE) 2 ]. (61) 

Here the key point to understanding a scalar ghost is that there is no more doipdoAE which 
is cancelled when integrating over B and (f> (see Eqs. (60) and (59)). Furthermore, the scalar 
mode E turns out to be non-dynamical since there is no time-derivative term of E, so we can 
integrate over the mode E. Consequently, becomes a Lagrangian for only a remaining 
scalar mode ip. A short calculation shows that 

= -^[ m *F xx + ^^-m 2n e- 2 ^ A }d ^ 
(A) 2 2n-l 6 

+ 2ne 2A ijAij-^--j^m 2n e- 2 ^ A ij 2 . (62) 

This Lagrangian gives us several important information on ip. First, as far as ghost is 
concerned, the scalar mode ip never be a ghost when I = ^-|^-[m 4 Fxx + ^pilf m 2ri e~ 2( "~ 2 ) A ] > 
0. Recall that Fxx > from one of the ghost condensation ansatzes (3), so / is definitely 
positive as long as n > \. Moreover, owing to the overall factor jj^ = 2{n — l) 2 {t — t ) 2 , I 

becomes divergent when t — > oo (we assume here that F xx is a positive constant as in the 
usual ghost condensation model), implying that the mode ip is certainly dynamical. Second, 
if we define the square of the mass of ip by 

Ml = -^——m 2n e - 2 ^ A , (63) 
^ 2n-l 3 ' K ' 
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it is also positive for n > ^ and approaches zero in the limit t — > oo like the tensor modes. 
The region of n > \ is consistent with n < 1 for the positive acceleration and ra > for 
the positivity of the square of the graviton mass. Third, as seen easily in (62), we have an 
unusal time-dependent dispersion relation, so the scalar mode ip in general breaks the Lorentz 
invariance. Finally, in the case of n — 0, this model reduces to the original ghost condensation 
one. 

Hence, in this section, we have explicitly shown that only the propagating modes around 
the cosmological backgrounds in this simple massive gravity are the massive tensor modes hjj 
and one scalar mode ip. As one peculiar feature, the massive scalar mode breaks the Lorentz 
symmetry manifestly. Even if we have limited ourselves to the simple massive gravity, we 
believe that this feature is also shared by the most general massive gravity models. 



5 The absence of non-unitary mode in the 't Hooft 
model with ghost condensation 

In this section, for completeness, we shall present a proof of the absence of the non-unitary 
mode in the 't Hooft model with ghost condensation which was considered in section 2. 

Now let us start by considering the fluctuations around the unitary gauge (8) for the four 
scalar fields 

(f) a = maTS* + ix a . (64) 

Then, X is expanded as 

X = -m 2 e- 2A + e~ 2A (-m 2 h 00 + 2md iv ) 

= X + e- 2A (-m 2 h 00 + 2md n ). (65) 

With Eq. (64), the 0*-equations of motion take the form 

- -&h - — Dtt 1 + 2A{h° l - — <9V) = 0. (66) 

2 Tfl Tft 



And the 0°-equation reads 



d [F xx (X)(h 00 - -d 7t )} = 0. (67) 
m 



At this stage, let us note that in the usual ghost condensation models which satisfy the 
ansatzes (3), one takes, for instance, the potential F(X) to be F(X) — c±(X — X) 2 + c 2 with 
some constants c\ > and c 2 . Then, it is natural to assume d F X x(X) = 0. With this 
assumption, the 0°-equation reduces to 

d (h 00 - -d n ) = 0. (68) 
m 
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This equation is easily solved to 



h 00 - -d n = C(x% (69) 
m 



where C(x l ) is an integration function depending on x l but not t. We then find that C(x l ) 
can be absorbed into the definition of ttq by redefining 

Ttl 

7i ^n + -tC(x l ). (70) 

Thus, we can take C(x l ) = 0. Then, n can be expressed in terms of h 0Q , so 7r is not an 
independent degree of freedom and can be neglected from the physical spectrum. 

The remaining equations of motion which we have to examine are Einstein's equations. 
At this point, for spatial diffeomorphisms, we take the gauge conditions 

7T* = 0. (71) 

With the help of the gauge conditions (71) and the 0°-equation (69) with C(x l ) = 0, Einstein's 
equations read 

l -[d p d,K + d p d v h? - □ V - d&h - v^(d P d a h^ - ah)] 

2 

772 

+ A[-2 Vlxl/ (d p h°» + ±doh) + 2r%(h)} 
m? 3 

= rrfSffiSij - —(3 - hu)r]^ - -m 2 h pu , (72) 

where Y° (h) denotes the Christoffel symbol in the linear approximation of h. 

In order to transform Einstein's equations (72) to the more tractable form, let us consider 
the following quantity 

P» = d v h v » - \d»h + 2Ah^. (73) 

With the ^-equations (66) and the gauge conditions (71), we have 

P i = 0. (74) 

Under the general coordinate transformations (57), the time-like component of P M transforms 

as 

5P° = DCo - 2i,9 Co - 2m 2 C . (75) 
Thus, using the remaining time-like diffeomorphism, we shall take a gauge 9 

P° = 0. (76) 



3 The same gauge condition was also taken in [10]. 
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Then, a little thought or calculation shows that Einstein's equations reduce to 



□ V - 2Ad Q V = 2m 2 V, (77) 

where we have defined /i MJ , = — ^rj^h. 

Notice that the right-hand side of 5P° in Eq. (75) has the same expression as the equations 
of motion for h^ v . Given that O( — 2Ado( — 2m 2 ( = 0, we obtain 5P° = 0. In other words, 
there remains a residual gauge symmetry associated with such the ( . The existence of the 
residual symmetry makes it possible to take a gauge 



h = 0, (76 



from which we have the equation 



hp/, = h^. (79) 
Finally, introducing H^ u = e^h^, it turns out that the whole equations read 

UH lxv --m 2 H lxv = 0, 

d u H^-AH 0fl = 0, 
H = 0, 

e- A H 00 - -d n = 0. (80) 
m 

These equations show that the graviton has mass \J~^m and the same 5 degrees of freedom 
as usual massive graviton modes. Moreover, there is no non-unitary mode since 7To mode is 
expressed in terms of H 00 . In this way, we find that the 't Hooft model with ghost condensation 
describes a physically plausible massive gravity model in the linearly expanding universe with 
zero acceleration. 



6 Conclusions and Discussion 

In this paper, we have shown that the 't Hooft model with ghost condensation is free of non- 
unitary scalar mode and is a massive gravity model in the linearly expanding cosmological 
universe with zero acceleration. This proof is rather general and simple. The reason is that in 
this model the acceleration is vanishing, from which many of equations take tractable expres- 
sions. Notice that the situation where there is no acceleration in the Friedmann-Robertson- 
Walker metric is similar to that where the equation of state is P = — \p with P and p being 
respectively the pressure and the matter density. This analogy might be useful for a better 
understanding of this solution. 

Furthermore, we have showed that a more general massive gravity model has an interesting 
class of classical solutions with the property of polynomially expanding cosmological universes 
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with non-zero acceleration. This class of solutions is classified by a constant n in the potential. 
Recall that this constant n is not a completely free parameter but receives some restriction 
from physical conditions. The first requirement that the model should describe a positive 
acceleration leads to n < 1. The second requirement that the massive graviton should not 
be a tachyon gives us a condition n > orn< \^ ( n !)• Moreover, the third 

requirement that the scalar mode should not be a ghost and/or a tachyon provides a final 
condition n > \. As a result, there exists the parameter region | < n < 1, which satisfies 
three requirements above at the same time. 

The existence of non-zero acceleration in the polynomially expanding cosmological uni- 
verses has made it difficult to prove that this massive gravity is free of the non-unitary mode 
and the Boulware-Deser instability. In order to clarify the physical propagating modes, we 
have used the (l+3)-parametrization of the metric fluctuations. Using this parametrization, 
it has been explicitly shown that the tensor modes become massive, the vector modes are non- 
dynamical, three of the scalar modes are also non-dynamical and one scalar becomes massive. 
This scalar is originally non-unitary mode, that is, a ghost, but becomes a normal particle 
because of the ghost condensation mechanism. However, the dispersion relation is not usual 
one, so this mode breaks the Lorentz invariance as expected from the ghost condensation 
scenario. 

In most of models which attempt to explain the late-time cosmic acceleration, the accel- 
eration is driven by some exotic matter with negative pressure called dark matter. On the 
other hand, in the models considered in this paper, the acceleration is driven by the massive 
graviton and an extra scalar which are originally part of components of the metric tensor 
h^y. This fact is in sharp contrast to the models which have been proposed so far. To our 
knowledge, no serious attempts have been made to study the late-time cosmic acceleration 
this way. 

As future's problems, we first wish to construct a Lorentz-invariant massive gravity model 
in a flat Minkowski space-time such that the model is free of the non-unitary mode since 
such a model describes a world of QCD. We also wish to examine various phenomenological 
aspects of the models that we have considered in this paper. It is known that the models lead 
to interesting phenomenology around a flat Minkowski background [14]- [16], so we think that 
the models around our cosmological backgrounds also give us new definite predictions for a 
scenario for inflation and density perturbations. 
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